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Abstract 

We report recent results of molecular dynamics simulations of fric- 
tional athermal particles at constant volume fraction and constant 
applied shear stress, focusing on a range of control parameters where 
the system first flows, but then jams after a time tj am . On decreas- 
ing the volume fraction, the mean jamming time diverges, while its 
sample fluctuations become so large that the jamming time probabil- 
ity distribution P{tj am ) becomes a power-law. We obtain an insight 
on the origin of this phenomenology focusing on the flowing regime, 
which is characterized by the presence of a clear correlation between 
the shear velocity and the mean number of contacts per particles Z, 
whereby small velocities occur when Z acquires higher values. 



1 Introduction 

The non-equilibrium transition from a fluid-like state to a disordered solid- 
like state, known as jamming transition, occurs in a wide variety of physical 
systems, such as colloidal suspensions and molecular fluids. Its widespread 
occurrence suggested to introduce a 'jamming phase diagram', with axes 
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temperature, density and shear stress [H [21 El E], able to describe the jam- 
ming transition of all of these systems. In this framework, jamming oc- 
curs at small temperature, high density, and small external shear stress. 
Frictionless athermal systems [1] can be described by the zero temperature 
plane of the jamming diagram: a transition line is expected to separate a 
jammed phase, which consists in non-ordinary solid Ej, and a flowing 
phase, where only steady states are observed [H |S]. Such a simple pic- 
ture becomes more involved in the presence of friction, which must be taken 
into account in order to properly describe granular (macroscopic) materials 

Ea nn nana mi sans- 

We have recently investigated the effect of friction on the jamming phe- 
nomenology performing Molecular Dynamics (MD) simulations of a frictional 
granular system at constant volume fraction and constant applied shear 
stress [151 US]- I n this set-up the jamming phenomenology is very rich, 
due to the occurrence of four possible phases, 'Flow', 'Flow & Jam', 'Slip & 
Jam' and 'Jam' (as found on increasing the density or decreasing the applied 
shear stress). In the 'Flow & Jam' and in the 'Slip & Jam' phase, the system 
reaches an apparently steady flowing state before jamming, or jams after a 
tiny displacement, respectively. In this paper, we present recent results re- 
garding the 'Flow & Jam' regime, where the system jams after flowing for 
a time tj am . We describe both dynamical properties of the system, such as 
the jamming time, as well as structural ones, such as the evolution of mean 
contact number. Connections and correlations between static and dynamics 
quantities suggest which is the physical mechanisms responsible of the 'Flow 
& Jam' phenomenology. 

2 The investigated system 

We perform MD simulations of the system considered in [T5l [T6] . Monodis- 
perse spherical grains of mass M and diameter D are enclosed in a box of 
dimension l x = l y = 16D, and l z = 8D. Periodic boundary conditions are 
used along x and y, while the size of the vertical dimension is fixed and chosen 
to be comparable to that of recent experiments [TTJ [TS1 US]- The upper and 
lower boundary surfaces of the box consist in rough "virtual" plates: each 
plate is made by a collection of particles that move as a rigid object. The 
bottom plate has an infinite mass, and is therefore fix, while the top one has 
a mass equal to the sum of the masses of its particles (roughly l x l y ). The top 
plate is subject to a shear stress a xz (a xz = a from now on). 

Grains interact via the standar linear spring-dashpot model. Two parti- 
cles % and j, in positions and r^, with linear velocities Vj and Vj, and angular 
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velocities cu, and ujj , interact if in contact, i.e., if the quantity £y = D — |r^-| 
is positive. 8y is called the penetration lenght, and = — is the dis- 
tance between the center of the particles i and j. The interaction force has a 
normal component F niJ and a tangential one F tij , the both having an elastic 
and a dissipative component: 

F nij = -kJijWij - 7„m e// v nij (1) 

Fty = -A; t u ti . - lt m eff v tij (2) 

where /c n and /c t are elastic moduli, while 7„ and 7f account for dis- 
sipative caharacter of the normal and tangential component rispectively. 
riij = Tij/lTijl \ nij = [(vi - Vj) • nijjriy, v ty = v y - v ni . and m e// is 
the effective mass. u t; ., set to zero at the beginning of a contact, measures 
the shear displacement during the lifetime of a contact. Its time evolution 
is fixed by Vt r and and u)j, as described in Ref. [20] • The presence of 
tangential forces implies the presence of torques, = — l/2rjj x F^. The 
shear displacement is set to zero both when a contact finish (Sij < 0), where 
fi is the coefficient of static friction. We use the value of the parameters 
of [20]: k n = 2 10 5 , k t /k n = 2/7, -y n = 50, ^t/ln = 0. Lenght, masses and 
times are expressed in units of D, m and \Jm/k n . 

The volume fraction (p represents the volume occupied by the grains di- 
vided by the volume of the container, i.e. <p = Nv /V +A(j), where V = l x l y l z 
is the volume of the system and vo = 1/QttD 3 is the volume occupied by a 
single grain and A0 = 0.021 is a corrective a term which takes into account 
the effect of the rough plates protruding into the system. By changing the 
number of particles we vary the volume fraction. 

The initial state is prepared setting to zero the friction coefficient |21j . 
randomly placing the particles and then inflating them until the desired vol- 
ume fraction is obtained; such a protocol is a short-cut of experimental pro- 
cedures with which it is possible to generate very dense disordered states of 
frictional systems, such as oscillations of high frequency and small ampli- 
tude [22]. After preparing the system we 'switch on' friction, and follow the 
time evolution of the system under the action of the constant shear stress . 
We investigate the evolution of the system focusing on a value of Coulomb 
friction coefficient equal to fi = 0.1 while (p and a are the variable control pa- 
rameter. All the data shown here refers to points (</?, a) falling in the 'Flow & 
Jam' region of the Jamming-phase-diagram [TB] . For each point we perform 
50-100 different runs, starting from different initial conditions. In the 'Flow 
& Jam' region, the steady state last a long time at small volume fractions, 
where the time tj am the system flows before jamming is large. This is the 
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Figure 1: (Color online) Velocity of the upper plate V s clS db function of time 
for three simulations at <fi ~ 0.627 and o = 2 10~ 3 . We marked the values of 
the jamming times, tj am , which differ by more than three decades. 




Figure 2: (Color online) Probability distribution P(tj am ) of the jamming 
time for a = 2 10 -3 and different values of the volume fraction, as indicated. 

region where the 'Flow & Jam' is more surprising, we have investigated in 
detail performing simulations lasting up to a time T = 5 10 4 . At higher vol- 
ume fraction tj am is small, and the identification of the steady state becomes 
more difficult. 
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3 Flow and Jam region 



3.1 Jamming times 

In the 'Flow & Jam' phase of sheared granular systems [16], the system 
first flows with a constant velocity as in a steady flowing phase, but then 
suddenly jams after a time tj am . Such a phenomenology is shown in Fig. [TJ 
where we report the velocity v s (t) of the upper plate as a function of time 
for three simulations at <p = 0.627 and a = 2 10 -3 . The average jamming 
time (tjam) depends on the volume fraction. It diverges as a power law on 
decreasing the volume fraction, at a critical volume fraction <f)ji(<J, 4>), and 
vanishes on increasing the volume fraction, at a critical value ^2(0", 4>)- The 
lines (f)ji(a,(f>) and </>j 2 (o", 0) define the boundary of the 'Flow & Jam' phase 
in the volume fraction, shear stress and friction jamming phase diagram for 
frictional particles [T6] . 

However, the jamming time is subject to large fluctuations. For instance, 
the three simulations shown in Fig. [TJjam at very different times, even though 
they differ only in the initial configuration. Such an observation suggested to 
study the sample fluctuations of the jamming time, which we show in Fig. [2] 
for a = 2 10~ 3 and three different values of the volume fraction. While at 
high volume fractions P(tj am ) is peaked, meaning that the jamming time 
is well defined, on decreasing the volume fraction the distribution moves to 
larger times, and at the same times changes shape, becoming well described 
by a power law. 

In order to understand the origin of the 'Flow & Jam' phenomenology 
and of the large fluctuations of the jamming time, which are not described 
by current rheological models such as those based on the inertial number [23] 
or on rate and state equations [211 [25], we have analyzed the evolution of the 
micro-structure of the system in the flowing regime and the micro-mechanics 
of the jamming process, as described below. 

3.2 Fluctuations of the micro-structure of the system 

Fig. E] illustrates the time evolution of the mean number of contacts per 
grain, Z(t), for the same simulations considered in Fig. [TJ At the beginning 
of the simulations, as a consequence of the considered preparation protocol, 
Z(t = 0) = 0. Z(t) rapidly increases as the system start flowing, and then 
fluctuates around a constant value in the following steady flowing phase. A 
comparison between Fig. [TJ and Fig. [3] suggests the presence of a correlation 
between the shear velocity v s (t) and the mean contact number Z(t), whereby 
large values of Z occurs when the shear velocity is small, and conversely. 
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Figure 3: (Color online) Mean number of contacts per grain Z as a function 
of time, for the simulations shown in FigCD The system jams when Z reaches 
a critical value Z c . 




Figure 4: (Color online) Parametric plot of the velocity of the top plate versus 
the mean number of contacts, obtained from the data shown in Fig.sdJand 
[3j The collapse indicates the existence of a correlation between Z and v s : 
the higher Z, the slower the system. 

Fig. H] shows a parametric plot of v s (t) versus Z(t). In such a representa- 
tion, the data from simulations characterized by very different values of the 
jamming time display a correlate behaviour. 

Note that the fluctuations of the shear velocity decreases as the mean 
number of contacts increases, suggesting the presence of a well defined critical 
mean number of contacts Z c at which jamming occurs, in agreement with the 
results of Fig. [3J We plot in Fig. \5\ the averaged shear velocity as a function 
of the mean number of contacts, (v s (Z)) ,The normal component is given by 
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Figure 5: (Color online) Dependence of the average shear velocity (v s ) as a 
function of the mean number of contacts per grain at a = 2 lO" 3 , for the 
indicated values of the volume fraction. 



(a) [.5 



0.5 






2 ior* 




2 L0~ 3 




2 icr 2 



KT 4.4 




Figure 6: (Color online) At = 0.629, a) average shear velocity (v s ) versus 
mean number of contacts per grain for the indicated values of the shear stress, 
and b) jamming critical mean contact number Z c as a function of the applied 
shear stress. 



for the indicated values of the volume fraction. The figure suggests that Z c is 
almost independent on the volume fraction, consistently with the existence 
of a volume fraction range where granular systems with equal mechanical 
properties can be prepared [16]. On the contrary, as shown in Fig.s[6k,b Z c 
increases with the applied shear stress. Considering that the shear modulus 
is expected to increase with the mean contact number [3] , this result suggests 
that the shear modulus of a jammed system does not simply depend on its 
volume fraction, but also on the applied stresses which caused jamming, or 
more generally on the preparation procedure. 
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4 Jamming mechanism 



Here we propose a qualitative mechanism to explain the origin of the 'Flow Sz 
Jam' phenomenology, based on the behavior of the shear velocity (v s (t)) and 
of the mean contact number Z(t), as well as on the dependence of tj am and 
Z c on the control parameters. The starting point is the well known dilatancy 
phenomenon observed in granular systems [23 ESI E2j , which is the tendency 
of flowing particulate systems to dilate. At constant pressure and constant 
shear strain rate a dilation is actually observed, the larger the greater the 
shear velocity [30]. At constant volume, which is the case considered here, 
dilation is obviously forbidden. This leads to a impeded dilatancy which may 
explain the observed phenomenology. 

While flowing, the system visits different microscopic configurations, each 
one having a typical mean number of contact Z. When Z is small, particles 
exert a small resistance to the applied stress, the shear rate increases and 
the system tries to dilate. This leads to a configuration with a larger mean 
number of contacts, exerting a larger resistance, which causes the system to 
decelerate. The existence of such a feedback mechanism is suggested by the 
correlations between Z and v s illustrated in Fig. HI The impeded dilatancy 
appears therefore responsible for the fluctuations of the mean number of 
contacts. The flowing system jams as a result of a large fluctuation of the 
mean number of contacts Z, which reaches the critical value Z c corresponding 
to configurations able to sustain the applied stress. How frequent are these 
fluctuations? We expect these fluctuations to be more rare when the volume 
fraction is small, simply because there are fewer configurations able to sustain 
the applied stress (i.e. with Z = Z c ): this explains why the jamming time 
increases as the volume fraction decreases. Also, one expects that when the 
volume fraction is smaller than a threshold value, there are no configurations 
with Z = Z c , which explains why the jamming time diverges decreasing the 
volume fraction. 

5 Conclusions 

In this manuscript, we focused on a region of the control parameters where 
frictional granular systems jam after flowing with a constant velocity, and 
described a possible mechanism able to explain the observed behavior, The 
mechanism is based on the notion of impeded dilatancy and on the presence 
of correlations between structure and dynamics of the system in the flowing 
regime. An open question ahead is the explanation of the power-law like 
distribution of the jamming times at small volume fractions, an indication of 
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the presence of correlations in the dynamics of sheared systems, which needs 
to be clarified. 
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